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1 Introduction 


Much work has been done recently to explore the SO q (3)-symmetric quantum 
mechanics developed in jlj and 0. In particular, a lot is known about the g- 
deformations of the harmonic oscillator. The other nontrivial soluble problem 
in ordinary quantum mechanics is the Coulomb problem, but for that one needs 
some notion of an inverse radius. Weich || considered a g-deformed Coulomb 
potential, defining | in a manner dependent upon a particular Hilbert space 
representation. This differs from the more standard ’’wave-function” type ap¬ 
proaches used in investigations of the oscillator (Refs. l-i for example). 

Here we approach this problem by defining as an actual element of the 
algebra of observables, thereby achieving representation-independence. Since 
X 2 = R 2 is already defined, we can then also define R as well as all its integral 
powers. A study of the action of momentum operators on powers of R then helps 
to bring out the interpretation of these operators as symmetric g-derivatives. 

Using this definition of -77, a self-adjoint symmetric g-deformation of the 
Coulomb Hamiltonian can be found which shares the n 2 -fold degeneracy of the 
undeformed Hamiltonian for the g-analog of bound states. As in 0, we obtain 
a Balmer-type spectrum for these states with 



where the symmetric g-analog of n is defined as 


n] q = 



( 2 ) 


In addition, we also obtain positive-energy wave functions and a candidate q- 
Coulomb S'-matrix. 

The paper is structured as follows. After a brief review of S'0,j(3)-symmetric 
quantum mechanics (Section 2), we set up the formalism for dealing with the 
g-Coulomb problem (Sections 3-5), which is then treated in Section 6. The 
detailed proofs of some statements made in the text are deferred to five Appen¬ 
dices. 


2 5'0 g (3)-symmetric quantum mechanics 

We build upon the algebra of observables as it is defined in Refs. [|| and 
The g-deformed metric and Levi-Civita tensors are defined as follows: 


T J = 7 ij = 


0 
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'0 0 -1 " 

£2 = Uj 2 = 0 — y/g 0 

_ 1 0 0 

0 0 o' 

0 0 75 

0 -y/q 0 J 

These give rise to an R-matrix: 

0 %j = qSlSj - e a «c fc , ° + (^ - l)7 ij '7fci 
which is a solution to the Yang-Baxter equation 

p/? pfr/c pac _ p?fc pid p/e 

^ab^cn^lm -^de^lf^rnni 

and has the inverse 

R - 1 Z = - Co ij Zkl a + (q- 

We use the metric and Levi-Civita tensors to define scalar and vector prod¬ 
ucts as for the undeformed tensors: A - B = y 1 ' 1 A, B 3 and [Ax B]^ = e*, 11 A t B 3 . 

4^ 3 is the S'0 9 (3)-covariant *-algebra defined by the generators X\, X 2 , X 3 
subject to the relations 

x X] k = 0 

and 

X* = '.'TV,. 

X 2 = X ■ X is then real and central in this algebra. In the <7=1 limit, X\ and 
X 3 correspond to -4^ {X ± Y) while X 2 is Z. The space of wave functions in 
harmonic-oscillator treatments of S'0 (} (3)-symmetric quantum mechanics is an 
appropriate subspace of A'^ 3 ■ 

One also considers an SO^Sj-covariant *-algebra 2>^ 3 of operators on A^ 3 . 

whose generators are the X,;, derivative operators 9,;, and a scaling operator /i. 
The Xi act on 44 3 by left multiplication, /1 is defined such that /i(l) = 1; and 

^q _ _ 

for all / € 4^ 3 , n{Xif) = qXip(f). 9* is defined such that 9,(1) = 0; and for 
all/Gx!f 3 , 

dM) = [ 7 a + ^x k di]f. 

The generators of V^ 3 then obey the relations: 

nXi = gXi^i, 
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/id; = 


-di, u, 

9 


diXj — 'jij H— R^jXkdi 


[5 x d] fe = 0. 


One also defines an inverse of /x: 

H- 1 = (j,[ 1 + g" 2 (l - q 2 )X ■ 3 + q~ 3 ( 1 - 9 ) 2 X 2 a 2 ]. 

In addition, there is a conjugate set of derivative operators in V ^ 3 , 

di = ii 2 [di + ( q~ 2 - 5 _ 1 )XjS 2 ]. 
di then satisfy the relations 

[9 x d\ k = 0 

and 

diX 0 = li3 + qR - 1 VXkdt 

The *-operation on /x and di is defined as fi* = q~ 3 fi~ x and (5,)* = —q 3 jddj. 

Neither triplet of derivative operators has a subalgebra isomorphic to ThXj, 
but a linear combination of the two does. This linear combination is then the 
triplet of (/-momentum operators 


p _ di + q 3 di 
* ~~ i( 1 + q~ 3 ) ' 

Then (Pi)* = 7 ^ Pj, [P x P\k = 0, and P 2 is a real scalar that commutes with 
the Pi. 

The Xi and Pj satisfy g-deformed versions of the Heisenberg relations: 
i(P a X b - qR- 1 %X c P d ) = /X - 1 (^TabW + ^-e ab 

and 

-i(X a P b - qR- 1 liP c X d ) = q 3 [i (^ ab W + ^-e ab m L m ^ , 
where W = /.x[l + g _2 (l — q)X ■ 9] is a real scalar, Lj = i/x[X x d\i , and 

p- = 9 -l + i. (3) 

9 
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The Li and W generate a g-deformed angular momentum algebra. Vectors Zi 
(such as the X ; t , dj, and Pk) satisfy the following relations with the L, and W, 
generalizing the role of the L,; as generators of rotations: 

LiZj = -a cd e dj e Z c L e + eij a Z a W, 


LZ=ZL= 0, 


WZj = KZjW - (K - l)ej rs Z r L s . 


In addition Z 2 commutes with the Li and W. 
following relations: 


L 2 = 


W 2 - 1 


K - 1 


The Li and W also satisfy the 


and 

[L x L] k = L k W = WL k . 


A ^3 modulo powers of X 2 can be shown to be a direct sum, indexed by 
nonnegative integers l, of irreducible representations of the angular momentum 
algebra. The Mr representation is then (21 + l)-dimensional, and IT is a Casimir 
operator with eigenvalue 


Wi 


qi+\ + q -l 

q + 1 


(4) 


For q = 1, the eigenvalue of L 2 in the /tli representation becomes the familiar 

1(1 + 1). 

As a last preliminary result, the momentum operators can be expressed in 
terms of W, /i, and X 2 : 


X 2 Pi 


iK (q 


1 


[-X.Wn- 1 

q 


WX i fi- 1 +q 2 X i Wn-qWXifi], (5) 


and 


X 2 P 2 = 


M 


1^2 1 


K 2 (4-j) 


q q 


( 6 ) 


These identities will be essential for calculations in the representation we will 
introduce later. We have included the numerical K factors on the left hand 
side of these identities, as we will define the hamiltonian scaled by a K 2 factor, 
which goes to one when q —> 1. 


3 Definition of 

A is already a well-defined concept in the space of undeformed, complex func¬ 
tions on R 3 . Its essential properties are that it is a real, scalar function and 
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that X 2 (±) 2 = (j^) 2 X 2 = 1. The simplest generalization of these properties is 
then to define ^ to be a real, scalar corepresentation of SO q { 3). is then 

the *-algebra generated by the X, and where the Xi obey the same relations 
as in A* 3 , commutes with the Xi, and 


X 2 {-) 2 = 1. 


(7) 


is the S'O g (3)-covariant *-algebra of operators on A ^, 3 , where we add 

■jj to the generators of V^ 3 . For this definition to be complete we must have a 

set of relations involving -b, and the derivatives. These must be consistent 
with the relations between -b and the Xi. Clearly since -b has dimensions of an 
inverse length, we should have 


»R =q 


-l 


ir 


Equation [t] implies di(j^) may be. 


d l X 2 (-) 2 = d z (l) = 0 
It follows from the algebra of T>^ 3 that 

d i X 2 = (q- 1 + l)X i + q 2 X 2 d i . 


Thus 


and therefore 


(q 1 +1 )Xijp + q 2 x 2 di(—j) — 0, 


R 2 ) ^ 


1 )— r T, 

) Ri 


( 8 ) 


Note that we are using the notation that if A £ 3 , and / £ Ar 3 then 

A(f) £ A^ 3 is the result of evaluating the effect of A on /, whereas Af £ 2>i 3 
is the product of A and / as operators. 

The simplest solution for di -g is 


ab^-'ba-^Tx,. 


Repeated application of Equation (j|) indeed gives 


-a 1 1 ~ 2 1 v 

“ K> x - 


(9) 

Similarly one finds that 

( 10 ) 


Equation (|j) must be checked for consistency with the algebra of A^ 3 before 
we can conclude that is a consistent operator algebra. In particular we must 


6 


show that di(j^f) = di(f^) and di(X 2 jp) = di(jpX 2 ) = 0. In addition if 
is truly a scalar, it should commute with the Li and W. The proofs that these 
conditions are satisfied are given in Appendix [A]. 

Having defined j we can now also define the g-deformed radius R = ^X 2 . 
This has the following commutation relation with <9^: 


d % R 


d ‘R X 2 


X 2 


q 1 + l)Xi + q 2 X 2 di] - q 2 -^Xi 

q~ 1 ^-X i + qRdi 
K 


Induction over positive and negative n gives 


d z R n = q-\n) q R n - 2 Xi + q n R n d,, 


where 


{n) q 


q n - 1 

q~ 1 


is the asymmetric g-analog. One can in fact develop a theory using d or d 
as momentum operators and rewrite g-deformed harmonic oscillator theories in 
the language of R and integral powers of R. However since d 2 and d are not 
self-adjoint, we will concentrate on the action of the Pi on elements of A^ 3 . 


4 Separation of Variables 

Using the formalism developed in the last section, we can now consider the q- 
analog of the separation of variables problem for the kinetic term of the hamil- 
tonian. To this end we first introduce g-analogs of spherical harmonics. These 
are defined, up to a—for us irrelevant—normalization as elements Y q l m of A^ :i 
which obey the following two conditions: 

(i) Left multiplication of Y q l m by L 2 or by L 2 yields Y q l m multipied by a real 
eigenvalue, or in other words, Y q l m is an “eigenfunction” of both L 2 and L 2 - 

(ii) All Y q commute with ^ l. 

For illustration, we give the explicit expressions of Y q l m for l = 0,1, 2 


II 

oo 

i 

ii 

i x 

y 1 _ 

I q o — 
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II 

i X3 

y 2 — 

1 q -2 — 


i m 

ii 

W XlX2 

V 2 — 

x q o 

[qX-i X: 

y 2 _ 

iql ~ 

— X 2 X 3 

y 2 — 

x q 2 — 

Y 


7f> 


Now, by multiplying equations (|j) and (j^) on the left by and jp respec¬ 
tively, we obtain 


KP, = 


i{q- ^)(q- !) R U 


and 


K 2 P 2 = 


1 1 


-XiW^- 1 - WXifjL- 1 + q 2 X, Wn - qWXifi 
(9 +l) 2 


- W 2 -qY~ 2 --( A ," 1 ) 2 

9 q 


(11) 

( 12 ) 


Since fx commutes with the Y q l m and W commutes with powers of R, if we ex¬ 
pand functions in terms of R n Y q l m , Equation ([l2|) is ready made for calculating 
their momentum squared. 


K 2 P 2 R n Y a L 


{q + l) \ 2 - q 2n+1 - 2 - g" 2 ”" 1 


1 1 

q 2l+l + q -2l-l _ q 2u+l _ q —2n—l , 

(9-i ) 2 W 


R n Y 1 

IX i q m 


In terms of the symmetric g-analog of n of equation (^), we can write this in 
the simplified form 

K 2 P 2 R n Y q l m = -[n + l + 1] Jn - l] q R n ~ 2 Y q l m . (13) 

This is a clear generalization of the result from ordinary real calculus that 

X[r n Y lm (d, $)} = (n + l + l)(n - l)r n ~ 2 Y lrn (6 , </>), 

and it is the main result of this section. 

Equation ( [Tl]) is not nearly as useful as its counterpart because the action 
of the X, on Yq is nontrivial. However for powers of R, one can obtain the 
simple result that 


iKPi(R n ) = [n\ q R n ~ 2 Xi = [n] q R 


yn-lXi 

R ' 


(14) 
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5 The free particle 


Let us consider a system with Hamiltonian 

H = K 2 P 2 


with the convenient normalization factor K 2 determined by equation (3). Then 
the Schrodinger equation for this system is the g-deformed Helmholtz equation 

K 2 P 2 ip = k 2 i/j. 

The solutions to this are of the form jt q u(kR)Y q and n\ q v(kR)Y q l m where j\ q u 
and ri[ q ]i are respectively the g-spherical Bessel and Neumann functions. 


3[q]l (*) 


OO 


E 

n=0 


(—l) w [2n + 2l] q \\ t 
[2?t] ^!! [2tz + 21 -\- 1]^! 


(15) 


i -1 

n w( x ) = E 

n =0 
oo 

+ (-i) m E- 

n=l 

where for nonnegative integers n 



and 



[2^ 2n]q! x 2n-l-i 


[2n] q U[2l - 2n] q \! X 


(—l) n [2n — 2l] q \\ 2 „_;_ x 
[2n] g !![2n - 2l\ q \ 

(16) 

nLiW? n >° 

1 n = 0 

(17) 

II 

5o 

JE 

►C5 

s s 

II V 

0 0 

(18) 


That Equations (El and © give rise to solutions is easily seen by applying 
Equation ©• This eliminates the zeroth element, and then reindexing gives 
the desired result as would be the case for the q = 1 differential equation. 

We can also obtain a g-deformed generalization of the Rayleigh formulas, 
which provide an alternative definition of these Bessel and Neumann functions. 
This will be a more convenient form for obtaining g-spherical Hankel functions, 


h [l]l( x ) = 3{q]i{x) + in [q] i{x) (19) 

and 

h {l] 2 ( X ) = J[q]l( X ) - in W 0 )> ( 20 ) 

corresponding to incoming and outgoing spherical waves. These results are 
discussed in Appendix [Bj. 

It is interesting to note that, because the Pi do not commute, it is not 
possible to have a plane wave with definite momentum as in ordinary quantum 
mechanics. The best we could do is specify the component of the momentum in 
the direction of propagation. Since, moreover, P 2 for example does not commute 
with L 2 , the problem of expanding even these quasiplane waves as a sum of 
R n Y q l m terms is nontrivial. 
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6 The g-Coulomb Problem 


In ordinary quantum mechanics, the Coulomb Hamiltonian is 

2m r 


or if we rescale this by 2m and incorporate the mass into a, 


H = p 2 


2a 

r 


There are several possible ways to g-deform this. We are interested in a self- 
adjoint Hamiltonian which preserves the properties of the ordinary Hamiltonian 
that make it amenable to finding eigenfunctions. That is to say we require the 
existence of a g-deformed Lenz vector which commutes with the g-deformed 
Hamiltonian so that there continue to be degeneracies between solutions with 
different angular momentum quantum numbers. 

Following ||, we define our g-deformed Coulomb Hamiltonian to be: 

H = K 2 P 2 ~aq(^ + p*^) 

which is clearly self-adjoint (for convenience the normalization of the kinetic 
term has again been chosen with a prefactor K 2 determined by equation (3)). 
Noting that //* = g _3 /x _1 , this can be written in the simpler form 

H = K 2 P 2 -a^{qp + q- 1 p,- 1 ). (21) 


This Hamiltonian also commutes with the Lenz vector 
. _ [W,iKP k ] aX k 

= - i - 

K-l R 
= iK(P k W - (P x L) k ) + 


The Lenz vector along with the angular momentum operators then generate the 
algebra given in Q. 

If we write the eigenvalues E of the Hamiltonian (|2|) in the form 


E = — 



( 22 ) 


then the 
tion) 


corresponding “eigenfunctions” are (see Appendix [c] for their 


'Ip'ylm — ^ ' dj, (y ) R. 

p =0 


, / aq 1 R\ p f —q l+1+p ~ 1 aR 




exp, 


W 5 


Y n 


Q m 


deriva- 

(23) 


with 


H p ( 7 ) = q pl +iP(P+ 1 ) 


(1 -q 2 ) p (g 2(i+1 ~ 7) ;g 2 ) p (g 4 ^ 1 ); q% 

(g 2 ; q 2 )p (g 2 ( 2 ' + 2 );g 2 ) P (g 2 ( i+ 1 );g 2 ) P ' 


(24) 
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Here 


( a;u) p = 


f nriotl-««”*) p =1,2,3... 

I 1 p = 0 


(25) 


is the g-deformed Pochhammer symbol j?|], with u = q 2 and q 4 in Eq. (^4j) , and 


exp 9 (x) 


OO 


E 

n =0 



(26) 


is the g-deformed exponential, where we have used the notation of Eq. (|l7|). g 7 
is obtained in terms of the energy E from the quadratic equation (^2|) . This has 
two solutions 


q 


- 7 ± 


V ± yV ~ E 


(27) 


with r) given by 



(28) 


(We assume g < 1 to insure convergence; had we instead set g > 1, we would 
have to everywhere change q —> g -1 ) 

Note that for a given energy, we actually only have one solution since if we 
write the solutions as power series of R, the difference equation admits only one 
solution that behaves as R l for small R. Whether we write it in terms of 7 + or 
7 _ simply gives us two expressions for the same result. 

In ordinary quantum mechanics, a decaying exponential multiplied by a 
polynomial in r is normalizable and gives rise to a bound state. Something very 
similar happens in the g / 1 limit. Consider the wave function (|23|), which 
is a sum of terms each of which is a nonnegative power of R multiplying a q- 
exponential of the form exp (—c q (p)R). As can be seen from (p3|), the c q (p) 
become independent of p as q —> 1, thus reproducing the usual result. Then we 
define a bound state in the g-deformed case to be a wave function of the type 
( p3| ) for which the sum over p truncates at the finite value p = n — l — 1. On 
account of the factor (g 2 d+ 1-7 ); g 2 ) p in the numerator o f (|24|) , such a truncation 
occurs if 7 equals a positive integer n > l. From Eq. (|22|), we then obtain the 
g-Balmer formula, and we find the corresponding wave functions given by ( p3| ) 
and (^dj) . The g-Balmer formula already appears in Ref. || , where, as noted in 
the introduction, the operator -^is treated differently. This difference is reflected 
in our wave functions. 

In principle these same wave functions (^3|) should also cover the continuum 
part of the spectrum, and one should be able to extract an S'-matrix from them. 
How to do this in a rigorous fashion remains to be seen. Here we consider the 
candidate 5-matrix suggested by equation (24) 


sl q) (E) 


1 q> r q2 (i + 1 - 7 _) 


00 


n 

n =0 


l _ g2(Z+l— 7 —bn) 

1 — g, 2 (Z+l— 7 ++n) ’ 


(29) 
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where the g-gamma function is defined as in [jjj: 

-A-- (») 

The ^-matrix appears to have all the right features: 

A. For q —> 1 this S'-matrix reproduces the familiar Coulomb S-matrix. In 

this limit the prefactor goes to one , as can be seen from Eq. and the q- 

gamma functions become precisely the ordinary gamma functions which appear 
in the ordinary Coulomb S'-matrix. 

B. For integer y + > l + 1, Si 9 \e) has a pole corresponding to a g-Balmer 
state. Both the the location and the residue of this pole differ from those of its 
ordinary (q = 1) Balmer limit. 

C. As can be seen from Eqs. ( p9| ) and (|?|) the branchpoint of S[ q \E), the 
scattering threshold, is now located E = rf, with 77 given in Eq. ( f28|) and not 
at E = 0 as in the ordinary case. This is the most dramatic departure from the 
ordinary case: the scattering region starts at E = rj 2 . For q = 1 this reduces to 
the expected threshold E = 0. 

A Properties of 

We have already satisfied di(X = 0. This definition of di also satisfies 
the other half of the second constraint. 

a ‘w = «- 2 w d ‘ ~ ^ + 




q 2 —^d i + q 2 X 2 d l 


[l)-q-\q^ + l)— Xl 


= 0 


Trivially <9,7 77 = di 

-^1 = di-g. Suppose that dif ^ = 

d t \f where / € 3 

= 

[7 ij+qH-^XbChlfl 


= 



= 

7 ijf^+qR-'VXkfq-'^ck- 

- Y ■) / 

" 

7 ijfft 

-q-^S^XkXtf + q- 1 ^ 

ia XkX t f-± 3 


= q-^ijf^ + ^R-'VXkdif 
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where we used the fact that e a kl XkXi = 0. At the same time 

= 9 _ 1 ^[ 7 « + ijX k di]f ~ q-^XiXjf 

= 

We must also show that <4 4 / 7i = Q i fLL ; n order to complete this inductive 
proof. However this is trivial having assumed that 9,;/A = Thus 

and "D^a are consistently defined. 

If -jj is truly a scalar, it should commute with the L.; and IT. This is indeed 
true: 


W = [l[l + q~ 2 {l-q)X -d] 

x-d = 7 ij Xidj^ 

= 4^3-4 

= i Ai [ 1 + (g- 2 - 9 - 1 )X.5] 


Li = liq x [X x d\i 

{Xxd] t = ei ik Xjd k 

= e4 Xj q—dk~q~-^Xk 

= ql-et^Xjdk 

Jtx 

= x ~ 5 \i 

Lij { = i-ur'[x x 4],;-^ 
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= -nq-'lXxdli 

1 , 

H Li 


B g-deformed Rayleigh formulas 


Let D be the symmetric g-derivative: 


= f(qx) - f(q 

{q — q~ x ) x 


(31) 


Acting on monomials, 

Dx n = [n] q x n ^ 1 (32) 

Let us define the following ^-generalizations of some common functions by 
replacing factorials with g-deformed factorials in their Taylor expansions: 


exp 9 (:r) = 

n —0 
oo 

(*) = H 


COS, 


n —0 


(-l) n x 2 
[2 n],! 


sin g (a:) = ^ 


(-D 


n^ 2 n+l 


n —0 


[2 n + l]qd 


Then exp g (z;r) = cos q (x) + i sin g (x) 
The g-deformed Rayleigh formulas, 

3{q]l{x) = i~ x ) 1 


-D 


sin, (a:) 


n[ q ]i(x) = (- x) l ( 


— C0Sg(a:) 


can then be proved by induction after noting that 

singer) 


j[q]o{ x ) = 


and 


n [q\0 X — 


COSg(x) 


The g-deformed spherical Hankel functions are 


h [l]i( x ) = J{q]i( x ) +in [q] i(x) 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 
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and h^y (x) is just the complex conjugate of h[H(x). Since is a 

linear operator, it follows that should satisfy 



(x) 


(-*)' 



^ -iexp q (ix) ^ 


One can show by induction that 


h[^(x) = qhW+P-^n+V)] 

71=0 


+ n\ q \ exp q (iq n x) 
[2 n] q \\\l — n\ q \ x n+1 


(41) 


satisfies this Rayleigh formula for all l. Thus these functions must equal the 
Hankel functions. 

The powers of q that appear both inside and outside the q-exponential arise 
because 


D(x k exp q (ax)) = [k] q x k 1 exp q (aqx) + aq k x k exp q (ax) 

= [— k\ q x k ~ x exp (a-a:) + aq k x k exp (ax). 

H q y 

This results from the g-deformed arithmetic in which 

[n + k] q = q n [k\ q + q~ k [n\ q = q~ n [k\ q + q k [n\ q (42) 

It appears to be a common trend of solutions to self-adjoint g-deformed Hamil¬ 
tonians that they can be expressed as series of the form )T) n A n exp q (q n x). 


C q -Coulomb Hamiltonian Eigenvalue Problem 

Let p = \J—E and [ 7 ]^ = (|. In order to obtain a difference equation for the A p , 
we need to express (H + /3 2 )ip as a series of the form J2 p B p exp q (q p x) where 
B p is a function of the A p . 

If we simply apply (H + /3 2 ) to tp using Equation [Td], we get 
(H + p 2 )R l+p exp q (-q s pR) Y q l m = 


g [n + p+2l+l] q [n + p\ q { _ qS p )nRn+p+l _^ ^ 

n=0 Wl- 




^n+Z+p—1 _|_ q—n—l—p—l 'j 


71=0 


[n] q \ 


(-q s /3) n R n+p+l - 1 Y q l m 


+(3 2 R l+p exp {-q s PR)Y q l m . 


If we decompose [n + p + 21 + l] q [n + p\ q , we get terms proportional to 
[n] 9 [n — l] g , q~ n [n] q , and q ~ 2n . By resumming these terms, we can rewrite 
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them as powers of R times exponentials. We wish to rewrite the entire equation 
in terms of functions of the form R l +P~ m exp q (—q s+m /3R) Y q l m where m is an 
integer and s is a function of p. Then we can obtain a difference equation for the 
coefficients of these functions. In the following, we assume that s = —l—l—p+ 7 . 
The only terms which need to be rewritten are then 


(K 2 P 2 ex Pq (—q ~ l ~ 1 ~ p+ ' 1 pR) Y q l m = 

[n + 2 l+p+l] q [n + p\ q 

( q l+i+P-lp) n R n + l +P~ 2 Y q l m 

^0 ^ 


-[7 ]g/?E 

n —0 


~n+Z+p+1 

1 (-g i+ 1 +P-T' i 9)" J R"+ i +P- 1 r g l m 

[ n W- 


= ~E ([n + 2 l+p+l] q [n+p] q -q 1+n+7 [n\ q [i] q ) 

n =0 


{-q l+ 1 +p-^py 


pn+Z+p —l 
1 Q m 


= - E {(q 2p+2l+2 + 1 - q 2l Mq[n - 1], + q~ 2n [p + 21 + 1 } q [p\ q 


n =0 


+ Q~ n [n] q (-q 1+1 bU P+2l+1 [p\ q + q 1+p [2l + 2 + p} q ) + q~ 2n [p + 21 + 1 ] q \p\ q } 

(-q l+ 1 +p -T/3) n 


pn+Z+p —l 

~ 11 1 q m 


J q- 


= + (g 2 *+ 2 + 2 P - q 2l+2+2p-2 1 _ q 4l+4p+4-2^ R l+p exp ^_ q -l- l-p+7^) ^ 

-/? (b +P+1 [+ - 9 2 p+ 3 i+ 1 ” 7 b ] 9 - ^+ 2 p+ 1 - 7 [ 2/ +p + 2 } q ) 

xR'+r - 1 exp q (—q~ l ~ p+ ' r pR) Y q l m 

-\p + 2i +1 } q \p\ q R l+p ~ 2 ex P g (-q~ l ~ p+1+ ' y PR) y q l m 

Thus 

(■H + P 2 )R 1+p ex Pq (-q- l - p - 1+ ^PR) Y q l m 
= ? 2,+2?+2 - 7 [2] r ,- P -i p + p +1 - q+'T 1 - q) 

xP 2 R 1 +p exp q (-q- l - 1 - p +'ypR)Y q l m 

+ (<? 2p+3/+1 - 7 b] 9 + q l+2p+1 ~^m +P + 2], - [2] g i+ P +i[qjg) 

x/JR'+P- 1 exp q (-q- l - p +^pR)Y q l m 
-\p + 2l + l} q [p} q R l+p ~ 2 exp q {-q- l ~ p+ ^ pR) Y q l m 
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If we sum over p from 0 to oo, we can obtain the desired difference equation. 
We normalize the wave functions so that Aq = 1. Then 


Ai = q l+1 (3 


[2] g i+i [l + 1 — 7] g 

[21 + 2] g 


and for p > 0 , 

{[p + 2Z + 3] g [p + 2] g } A p+2 
= « 2 ' +2p+2_7 [2] ? i+ P +i [Z + p + 1 - l\ q (q~ l - g)/? 2 A> 

+ (g 2 P+ 3 '+ 3 -7[p + 1], + g«+2p+3-7 [2i + p + 3]g _ [2] rl _ p _ 2 [ 7 ],) /M p+1 

Equation ^4] is then a solution to this difference equation. Thus the theorem is 
proved. 
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